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Introduction to new 
Hierarchic theory of condensed matter 
Qhttp: / /arXiv.org/abs/physics/00030044|) 

A basically new hierarchic quantitative theory, general for solids 
and liquids, has been developed. 

It is assumed, that unharmonic oscillations of particles in any con- 
densed matter lead to emergence of three-dimensional (3D) superpo- 
sition of standing de Broglie waves of molecules, electromagnetic and 
acoustic waves. Consequently, any condensed matter could be con- 
sidered as a gas of 3D standing waves of corresponding nature. Our 
approach unifies and develops strongly the Einstein's and Debye's 
models. 

Collective excitations, like 3D standing de Broglie waves of molecules, 
representing at certain conditions the mesoscopic molecular Bose con- 
densate, were analyzed, as a background of hierarchic model of con- 
densed matter. 

The most probable de Broglie wave (wave B) length is deter- 
mined by the ratio of Plank constant to the most probable impulse 
of molecules, or by ratio of its most probable phase velocity to fre- 
quency. The waves B are related to molecular translations (tr) and 
librations (lb). 

As the quantum dynamics of condensed matter does not follow in general 
case the classical Maxwell-Boltzmann distribution, the real most probable de 
Broglie wave length can exceed the classical thermal de Broglie wave length and 
the distance between centers of molecules many times. 

This makes possible the atomic and molecular Bose condensation in solids 
and liquids at temperatures, below boiling point. It is one of the most important 
results of new theory, which we have confirmed by computer simulations on 
examples of water and ice. 

Four strongly interrelated new types of quasiparticles (collective excita- 
tions) were introduced in our hierarchic model: 
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1. Effectons (tr and lb), existing in "acoustic" (a) and "optic" (b) states 
represent the coherent clusters in general case; 

2. Convertons, corresponding to interconversions between tr and lb types of 
the effectons (flickering clusters); 

3. Transitons are the intermediate [a ^ b] transition states of the tr and lb 
effectons; 

4. Deformons are the 3D superposition of IR electromagnetic or acoustic 
waves, activated by transitons and convertons. 

Primary effectons (tr and lb) arc formed by 3D superposition of the 
most probable standing de Broglie waves of the oscillating ions, atoms or 
molecules. The volume of effectons (tr and lb) may contain from less than one, 
to tens and even thousands of molecules. The first condition means validity 
of classical approximation in description of the subsystems of the effectons. 
The second one points to quantum properties of coherent clusters due to 
molecular Bose condensation. 

The liquids are semiclassical systems because their primary (tr) effectons 
contain less than one molecule and primary (lb) effectons - more than one 
molecule. The solids are quantum systems totally because both kind of their 
primary effectons (tr and lb) are molecular Bose condensates. These conse- 
quences of our theory are confirmed by computer calculations. 

The 1st order [gas — ► liquid] transition is accompanied by strong decreasing 
of rotational (librational) degrees of freedom due to emergence of primary (lb) 
effectons and [liquid — > solid] transition - by decreasing of translational degrees 
of freedom due to Bose-condensation of primary (tr) effectons. 

In the general case the effecton can be approximated by par- 
allelepiped with edges corresponding to de Broglie waves length in 
three selected directions (1, 2, 3), related to the symmetry of the 
molecular dynamics. In the case of isotropic molecular motion the 
effectons' shape may be approximated by cube. 

The edge-length of primary effectons (tr and lb) can be considered 
as the "parameter of order". 

The in-phase oscillations of molecules in the effectons correspond to the 
effecton's (a) - acoustic state and the counterphase oscillations correspond to 
their (b) - optic state. States (a) and (b) of the effectons differ in potential 
energy only, however, their kinetic energies, impulses and spatial dimensions - 
are the same. The 6-state of the effectons has a common feature with Frolich's 
polar mode. 

The (a — > b) or (b — ► a) transition states of the primary effectons 
(tr and lb), defined as primary transitons, are accompanied by a 
change in molecule polarizability and dipole moment without density 
fluctuations. At this case they lead to absorption or radiation of IR 
photons, respectively. 

Superposition (interception) of three internal standing IR pho- 
tons of different directions (1,2,3) - forms primary electromagnetic 
deformons (tr and lb). 
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On the other hand, the [lb^ tr] convertons and secondary transitons are 
accompanied by the density fluctuations, leading to absorption or radiation of 
phonons. 

Superposition resulting from interception of standing phonons in three direc- 
tions (1,2,3), forms secondary acoustic deformons (tr and lb). 

Correlated collective excitations of primary and secondary effectons and 
deformons (tr and lb), localized in the volume of primary tr and lb electromag- 
netic deformons, lead to origination of macroefFectons, macrotransitons 
and macrodeformons (tr and lb respectively). 

Correlated simultaneous excitations of tr and lb macroeffectons in the vol- 
ume of superimposed tr and lb electromagnetic deformons lead to origination 
of supereffectons. 

In turn, the coherent excitation of both: tr and lb macrodeformons and 
macroconvertons in the same volume means creation of superdeformons. Su- 
perdeformons are the biggest (cavitational) fluctuations, leading to microbub- 
bles in liquids and to local defects in solids. 

Total number of quasiparticles of condensed matter equal to 4!=24, 
reflects all of possible combinations of the four basic ones [1-4], intro- 
duced above. This set of collective excitations in the form of "gas" of 
3D standing waves of three types: de Broglie, acoustic and electro- 
magnetic - is shown to be able to explain virtually all the properties 
of all condensed matter. 

The important positive feature of our hierarchic model of matter is that it 
does not need the semi-empiric intermolecular potentials for calculations, which 
are unavoidable in existing theories of many body systems. The potential energy 
of intermolecular interaction is involved indirectly in dimensions and stability 
of guasiparticles, introduced in our model. 

The main formulae of theory are the same for liquids and solids 
and include following experimental parameters, which take into ac- 
count their different properties: 

[1]- Positions of (tr) and (lb) bands in oscillatory spectra; 

[2]- Sound velocity; 

[3]- Density; 

[4]- Refraction index (extrapolated to the infinitive wave length of 
photon). 

The knowledge of these four basic parameters at the same temperature and 
pressure makes it possible using our computer program, to evaluate more than 
300 important characteristics of any condensed matter. Among them are such 
as: total internal energy, kinetic and potential energies, heat-capacity and ther- 
mal conductivity, surface tension, vapor pressure, viscosity, coefficient of self- 
diffusion, osmotic pressure, solvent activity, etc. Most of calculated parameters 
are hidden, i.e. inaccessible to direct experimental measurement. 

The new interpretation and evaluation of Brillouin light scattering and 
Mossbauer effect parameters may also be done on the basis of hierarchic the- 
ory. Mesoscopic scenarios of turbulence, superconductivity and superfluity are 
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elaborated. 

Some original aspects of water in organization and large-scale dynamics of 
biosystems - such as proteins, DNA, microtubules, membranes and regulative 
role of water in cytoplasm, cancer development, quantum neurodynamics, etc. 
have been analyzed in the framework of Hierarchic theory. 

Computerized verification of our Hierarchic concept of matter on 
examples of water and ice is performed, using special computer pro- 
gram: Comprehensive Analyzer of Matter Properties (CAMP, copy- 
right, 1997, Kaivarainen) . The new optoacoustic device, based on this 
program, with possibilities much wider, than that of IR, Raman and 



Brillouin spectrometers, has been proposed ( http:/ /www. karelia.ru/~alexk 
[CAMP]). 

This is the first theory able to predict all known experimental 
temperature anomalies for water and ice. The conformity between 
theory and experiment is very good even without any adjustable pa- 
rameters. 

The hierarchic concept creates a bridge between micro- and macro- 
phenomena, dynamics and thermodynamics, liquids and solids in 
terms of quantum physics. 



1. Role of inter-domain water clusters in large-scale dynamics of 

proteins 

The functioning of proteins, namely antibodies, enzymes, is caused 
by the physicochemical properties, geometry and dynamics of their 
active sites. The mobility of an active site is related to the dynamics 
of the residual part of a protein molecule, its hydration shell and the 
properties of a free solvent. 

The dynamic model of a protein proposed in 1975 and supported 
nowadays with numerous data (Kaivarainen, 1985,19896), is based on 
the following statements: 

1. A protein molecule contains one or more cavities or clefts ca- 
pable to large scale fluctuations - pulsations between two states: 
"closed" (A) and "open" (B) with lesser and bigger accessibility to 
water. 

The frequency of pulsations (va^b)'- 



ioV 1 < i/ te B < 10V 1 

depends on the structure of protein, its ligand state, temperature and 
solvent viscosity. Transitions between A and B states are the result of 
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the relative displacements of protein domains and subunits forming 
the cavities; 

2. The water, interacting with protein, consists of two main frac- 
tions. 

The 1st major fraction, which solvates the outer surface regions of 
protein has less apparent cooperative properties than the 2nd minor 
fraction confined to "open" cavities. Water molecules, interacting 
with the cavity in the "open" (B)-state form a cooperative cluster, 
whose lifetime (> 10 10 s x ). Properties of clusters are determined by 
the geometry, mobility and polarity of the cavity, as well as by tem- 
perature and pressure. 

It is seen from X-ray structural data that the protein cavities: active sites 
(AS), other interdomain clefts, the space between subunits of oligomeric pro- 
teins, have a high nonpolar residues content. In contrast to the small intra 
domain holes isolated from the outer medium, which sometimes contain several 
H 2 molecules, the interdomain and intersubunit cavities can contain several 
dozens of molecules (Fig. 1), exchanging with bulk water. 

The development of the above dynamic model has lead us to the 
following classification of dynamics in the native globular proteins. 

1. Small-scale (SS) dynamics: 

low amplitude (< 1 A) thermal fluctuations of atoms, aminoacids residues, 
and displacements of a-helixes and /3-structures within domains and subunits, 
at which the effective Stokes radius of domains does not change. This 
type of motion, related to domain stability, can differ in the content of A and 
B conformers (Fig. 1, dashed line). The range of characteristic times at SS 
dynamics is (10~ 4 — 10~ n )s, determined by activation energies of corresponding 
transitions. 

2. Large-scale (LS) dynamics: 

is subdivided into LS-pulsations and LS-librations with a character 
of limited diffusion of domains and subunits of proteins: 

LS- pulsations are represented by relative translational-rotational displace- 
ments of domains and subunits at distances > 3A Thus, the cavities, which arc 
formed by domains, fluctuate between states with less (A) and more (B) water- 
accessibility. The life-times of these states depending on protein structure and 
external conditions are in the limits of (l(T 4 - 1CT 7 ) s. 

In accordance to our model, one of contributions to this time is determined 
by frequency of excitations of [Ib/tr] macroconvertons. The frequency of 
macroconvertons excitation at normal conditions is about 10 7 (1/s). 

The pulsation frequency of big multi-subunit oligomeric proteins of about 
10 4 (1/s) could be related to stronger fluctuations of water cluster in their central 
cavity like macrodeformons or even superdeformons (Fig. 3c, d). 

The life-times of (A) and (B) conformer markedly exceeds the time of tran- 
sitions between them ~ (10~ 9 — 10 -11 ) s. 

The (A B) pulsations of various cavities in proteins could be correlated. 
The corresponding A and B conformers have different Stokes radii and effective 
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volume. 

The geometrical deformation of the inter-subunits large central cavity of 
oligomeric proteins and the destabilization of the water cluster located in it 
lead to relaxational change of (^4 B) equilibrium constant: 

r , / Ga - Gb 

Ka^-b = exp - 



RT 

The dashed line means that the stability and the small-scale dynamics of do- 
mains and subunits in the content of A and B conformers can differ from 
each other. The [A ^ B] pulsations are accompanied by reversible sorption- 
desorption of (20 — 50) H 2 molecules from the cavities. 

Structural domains are space-separated formations with a mass of (1 — 2) • 
10 3 D. Protein subunits (MM > 2 • 10 3 D), as a rule, consist of 2 or more 
domains. The domains can consist only of a or only of /3-structure or have no 
like secondary structure at all (Schulz, Schirmer, 1979). 

The shift of A B equilibrium of central cavity of oligomeric proteins 
determines their cooperative properties during consecutive ligand binding in 
the active sites. Signal transmission from the active sites to the remote regions 
of macromolecules is also dependent on (A B) equilibrium. 




Fig. 1. Examples of large-scale (LS) protein dynamics: A B 
pulsations and librations with correlation times (r^ < (Kaivaralnen, 
1985, 1989): 

a) mobility of domains connected by flexible hinge or contact 
region, like in the light chains of immunoglobulins; 

b) mobility of domains that form the active sites of proteins, like 
in hexokinase, papain, pepsin, lysozyme etc. due to flexibility of 
contacts; 



7 



c) mobility of subunits forming the oligomcric proteins like hemoglobin. 
Besides transitions of the active sites of each subunit, the (^4 4=> B) 
pulsations with frequencies of (10 4 — 10 6 ) s _1 are pertinent to the 
common central cavity. 

b) Vibrations represent the relative rotational - translational mo- 
tions of domains and subunits in composition of A and B conformers 
with correlation times tm — (1 — 5) ■ 10~ 8 s. 

LS - librations of domains are accompanied by "flickering" of water 
cluster in the open cavity between domains or subunits. The process 
of water cluster "flickering", i.e. [dissociation ^association] is close 
to the reversible first-order phase transition, when: 

AGh 2 o — &-Hh 2 o — TASh 2 o ~ 

Such type of transitions in water-macromolecular systems could 
be responsible for so called "enthalpy-entropy compensation effects" 
(Lumry and Biltonen, 1969). 

The "flickering clusters" means excitation of [Ib/tr] conversions be- 
tween librational and translational primary water effectons, accom- 
panied by [association/dissociation] of coherent water cluster (see dif- 
ference in dimensions of lb and tr effectons on Fig. 18a, b of [1]). 

The water cluster (primary lb effecton) association and dissocia- 
tion in protein cavities in terms of mesoscopic model represent the 
(etc) - convertons or (be) - convertons. These excitations stimulate the 
LS- librations of domains in composition of B-conformer. The fre- 
quencies of (ac) and (be) convertons, has the order of about 10 8 c _1 . 
This value coincides well with experimental characteristic times for 
protein domains librations. 

The (ac) and (be) convertons represent transitions between similar 
states of primary librational and translational effectons: [an, ^ a tr ] and 
[bu, ^ hr] (see Introduction to[l, 2]). 

For the other hand, the Macroconvertons, representing simultane- 
ous excitation of (ac + be) convertons, are responsible for [B ^ ^4] large- 
scale pulsations of proteins. 

The librational mobility of domains and subunits is revealed by the fact that 
the experimental value of tm is less than the theoretical one (t 1 m ) calculated on 
the Stokes-Einstein formula: 

rh = (V/k) • V /T 

This formula is based on the assumption that the whole protein can be 
approximated by a rigid sphere. It means, that the large-scale dynamics can be 
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characterized by the " flexibility factor" , in the absence of aggregation equal to 
ratio: 

fl = (tm/tI) < 1 

Antonchenko (1986) has demonstrated, using the Monte-Carlo method for sim- 
ulations, that the disjoining pressure of a liquid in the pores onto the walls 
changes periodically depending on the distance (L) between the limiting sur- 
faces. If the water molecules are approximated by rigid globes, then the maxima 
of the wedging pressure lie on the values of distance L: 9.8; 7; and 3.3A. It points, 
that small changes in the geometry of cavities can lead to significant changes in 
their A B equilibrium constant (Ka&b)- 

According to our model the large-scale transition of the protein cavity from 
the "open" B-statc to the "closed" A-state consists of the following stages: 

1. Small reorientation (libration) of domains or subunits, which form an 
"open" cavity (B-statc). This process is induced by (ac) or (be) convertons of 
water librational effecton, localized in cavity (flickering of water cluster); 

2. Cavitational fluctuation of water cluster, containing (20—50) H 2 molecules 
and the destabilization of the B-state of cavity as a result of [lb^tr] macro- 
converton excitation; 

3. Collapsing of a cavity during the time about 10~ 10 s, dependent on pre- 
vious stage and concomitant rapid structural change in the hinge region of 
interdomain and intersubunit contacts: [B — ► A] transition. 

The b — > a transition of one of the protein cavities can be followed by similar 
or the opposite A — > B transition of the other cavity in the macromolecule. 

It should be noted that the collapsing time of a cavitation bubble with the 
radius: r ~ (10 — 15) A in bulk water and collapsing time of interdomain cavity 
are of the same order: At ~10~ 10 s under normal conditions (Shutilov, 1980). 

If configurational changes of macromolecules at B — > A and A — > B transi- 
tions are sufficiently quick and occur as a jumps of the effective volume, they 
accompanied by appearance of the shock acoustic waves in the bulk medium. 

When the cavitational fluctuation of water in the "open" cavity does not 
occur, then (b — > a) or (B — > A) transitions are slower processes, determined 
by continuous diffusion of domains and subunits. This happens when [lb — > tr] 
macroconvertons are not excited. 

In their review, Karplus and McCammon (1986) analyzed data on alco- 
holdehydrogenase, myoglobin and ribonuclease, which have been obtained using 
molecular dynamics approach. It has been shown that large-scale reorienta- 
tion of domains occur together with their deformation and motions of a and (3 
structures. 

It has been shown also (Karplus and McCammon, 1986) that activation 
free energies, necessary for [A B] transitions and the reorganization of hinge 
region between domains, do not exceed (3-4) kcal/mole. Such low values were 
obtained for proteins with even rather dense interdomain region, as seen from 
X-ray data. The authors explain such low values of activation energy by the 
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fact that the displacement of atoms, necessary for such transition, does not 
exceed 0.5 A, i.e. they are comparable with the usual amplitudes of atomic 
oscillation at temperatures 20 — 30° C. It means that they occur very quickly 
within times of 10~ 12 s, i.e. much less than the times of [A — > B] or [B — > A] 
domain displacements (10~ 9 — 10~ 10 s). Therefore, the high frequency small- 
scale dynamics of hinge is responsible for the quick adaptation of hinge geometry 
to the changing distance between the domains and for decreasing the total 
activation energy of [A B] pulsations of proteins,. 

Recent calculations by means of molecular dynamics reveal that the oscilla- 
tions in proteins are harmonic at the low temperature (T< 220K) only. At the 
physiological temperatures the oscillations are strongly unharmonic, 
collective, global and their amplitude increases with hydration (Stein- 
back et al., 1996). Water is a "catalyzer" of protein unharmonic 
dynamics. 

It is obvious, that both small-scale (SS) and large-scale (LS) dy- 
namics, introduced in our model, are necessary for protein function. 
To characterize quantitatively the LS dynamics of proteins, we pro- 
posed the unified Stokes-Einstein and Eyring-Polany equation. 



2. Description of large-scale dynamics of proteins, based on 
generalized Stokes-Einstein and Eyring-Polany equation 

In the case of the continuous Brownian diffusion of a particle, the rate 
constant of diffusion is determined by the Stokes-Einstein law: 



t Vr\ 

where: r is correlation time, i.e. the time, necessary for rotation of a particle 
by the mean angle determined as <p rts 0.5 of the turn or the characteristic time 
for the translational movement of a particle with the radius (a) on the distance 
{h x f/ 2 ~ 0.6a (Einstein, 1965); 

V = 47ra 3 /3 is the volume of the spherical particle; ks is the Boltzmann 
constant, T and r] are the absolute temperature and bulk viscosity of the solvent. 

On the other hand, the rate constant of [A — > B] reaction for a molecule in 
gas phase, which is related to passing through the activation barrier G A ~* B , is 
described with the Eyring-Polany equation: 



To describe the large-scale dynamics of macromolcculcs in solution related to 
fluctuations of domains and subunits (librations and pulsations), an equation is 
needed which takes into account the diffusion and activation processes simulta- 
neously. 
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The rate constant for the rotational- translational diffusion of the parti- 
cle (k c ) forming a macromolecule (continuous LS-dynamics) is expressed with 
the generalized Stokes-Einstein and Eyring-Polany equation (Kaivarainen and 
Goryunov, 1987): 

k B T ( G st \ _! 

where: V is the effective volume of domain or subunit, which are capable to 
the Brownian mobility independently from the rest part of the macromolecule, 
with the probability: 



,4, 

where: G s t is the activation energy of structural change in the contact (hinge) 
region of a macromolecule, necessary for independent mobility of domain or 
subunit; r c is the effective correlation time for the continuous diffusion of this 
relatively independent particle. 

The effective volume V can be changed under the influence of 
temperature, perturbants and ligands. 

The generalized Stokes-Einstein and Eyring-Polany equation (3) 
is applicable also to describing the diffusion of the whole (integer) 
particle, dependent on the surrounding medium fluctuations with ac- 
tivation energy (G a ). The ligand diffusion in the active site cavity of 
proteins is such a type of processes. 

To describe noncontinuous process, the formula for rate constant (fcj ump ) 
of the jump-like translations of particle, related to emergency of cavitational 
fluctuations (holes) near the particle was proposed (Kaivarainen and Goryunov, 
1987): 



r jump 



RT 



where: 



W = aS + n s (/i ou t - Min) (6) 

is the work of cavitation fluctuation with the cavity surface S, at which n s 
molecules of the solvent (water) change its effective chemical potential from 

to fl out - 

The dimensions of cavity fluctuation near particle must be comparable to 
corresponding particles. 
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In a homogeneous phase (i.e. pure water) under equilibrium conditions we 
have: /j, ln = /i ut- With an increase of particle sizes, surface of cavitational 
fluctuation (S) and its work (W), the corresponding probability of cavitation 
fluctuations: 

Pjump = exp(-W/RT) 

will fall. 

The notion of the surface energy (a) retains its meaning even at very small 
"holes" because of its molecular nature (see Section 11.4 of [1] and [4]). 

Tj™™ p in eq. (5) is the minimal possible jump-time of a particle with mass 
(m) over the distance A with the mean velocity: 

« max - (2fcT/m) 1 / 2 (7) 
Hence, we derive for the maximal jump-rate at W=0: 



Vmax 1 / 2kT^ 



jump 



7,max _ _ • max _ — ,„\ 

jump " r-l"„ A ~ A \ to ' 



In the case of hinged domains, forming macromolecules their relative A f± B 
displacements (pulsations) are related not only to possible holes forming in 
the interdomain (intersubunit) cavities or near their outer surfaces, but to the 
structural change of hinge regions as well. 

If the activation energy of necessary structure changes is equal to G^^ B , 
then eq. (5), with regard for (8), is transformed into 



Jump A V m J * \ RT 

where: Wa,b is the work required for cavitational fluctuations of water; this 
work can be different in two directions: (Wb) is necessary for nonmonotonic 
B — > A transition; (Wa) is necessary for jump- way A — > B transition. 

Under certain conditions A ^ B transitions between protein conformers 
(LS- pulsations) can be realized owing to the jump- way and continuous types 
of relative diffusion of domains or subunits as two stage reaction. In this case, 
the resulting rate constant of the process will be expressed through (9) and (3) 
as: 



h A^B _ h A^B _ h AoB _ 1 (2kT\ 1 / 2 pYn ( Wa,b+G^ b \ 
«Tes — K-jump ~ K-c A V m / ^ ^ RT ) ^ 

(10) 

-f- xv exp ^ RT j 
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The interaction between two domains in A-conformer can be described 
using microscopic Hamaker - de Bour theory. One of the contributions into 
Gf^ B is the energy of dispersion interactions between domains of the radius 
(a) (Kaivaralnen, 1989b): 

[Gst ~U H ~ -A*a/12H] AiB (11) 

where 

A* « {A l J 2 - Al' 2 ) 2 « ^Khv s [a s N s - a c N c } 2 (12) 

is the complex Hamaker constant; H is the slit thickness between domains in A- 
state; A c and A s arc simple Hamaker constants, characterizing the properties 
of water in the A-state of the cavity and in the bulk solvent, correspondingly. 
They depend on the concentration of water molecules (N c ~ N s ) and their 
polarizability (a c ^ a s ): 

A c = '^nhvQalN 2 and A s = ^nhvQa 2 s N 2 ; 

where: hv^ w hv$ are the ionization potentials of H2O molecules in a cavity 
and in a free solvent. 

In the "closed" A-state of a cavity the water layer between domains has a 
more compact packing as compared with the ice-like structure of a water cluster 
in the B-state of a cavity, or with a free solvent. As far Ha < Hb the dispersion 
interaction (11) between domains in A- state of cavity is stronger, than that in 
B-state: Ufi>Ug. 

Disjoining pressure of water in the cavities 

n = -A*/6nH 3 (12a) 

decreases with the increase of the complex Hamaker constant (A* ) that corre- 
sponds to the increase of the attraction energy (Uh) between domains. 

Cooperative properties of clusters in open (B)-states of the cavities are more 
pronounced as compared to that in bulk water. That results in the greater 
changes of a c N c than that of a s N s induced by temperature. The elevation of 
temperature decreasing the dimensions of intcrdomain water clusters leads to 
the strengthening of intcrdomain interaction, while the lowering temperature 
leads to opposite effect. 

We can judge about the changes of a s N s in the experiment on measuring 
the solvent refraction index, as far from our theory of refraction index (eq. 8.14 
of [1] or paper [3]): 

(n 2 - l)/n 2 s = \i:a s N s or: a s N s = jn ■ (12b) 
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In the closed cavities the effect of temperature on water properties is lower 
as compared to that in bulk water. It follows that thermoinduced non- 
monotonic transition in the solvent refraction index must be accom- 
panied by in-phase nonmonotonic changes of the [A <^> B] equilibrium 
constant (Ka^b)- As far (A) and (B) conformers usually have different sta- 
bility and flexibility the changes of Ka^b will be manifested in the changes 
of protein large-scale and small-scale dynamics. It has been shown before that 
viscosity itself has nonmonotonic temperature dependence due to the nonmono- 
tonic dependence of n 2 (t) (eg. 11. 44, 11.45 and 11.48 of [1] or paper [4]). 

Thus, thermoinduced non-denaturational transitions of macromolecules 
and supramolecular systems located in the aqueous environment are 
caused by nonmonotonic changes in solvent properties, including its 
refraction index. 

The influence of D2O and other pcrturbants on protein dynamics is explained 
in a similar way. The effect of deuterium oxide (D2O) is a result of substitution 
of H 2 from protein cavities and corresponding change of complex Hamaker 
constant (12). 

Generalized equation (3) is applicable not only for evaluating the frequency 
of macromolecules transition between A and B conformers but also for the 
frequency of the dumped librations of domains and subunits within A and B 
conformers. Judging by various data (Kaivarainen, 1985, 19896), the interval of 
A ^ B pulsation frequency is: 



V A ^B = -^T- « = ( 104 - 107 ) S_1 ( 13 ) 

M + tB 

where: tA and ts are the lifetimes of A and B conformers. 

The corresponding interval of the total activation energy of the jump-way 
A ^=> B pulsations can be evaluated from the eq. (9). We assume for this end 
that the pre-exponential multiplier is about 10 10 s _1 as a frequency of cavita- 

o 

tional fluctuations in water with the radius "(10-15) A. 
Taking a logarithm of (9) we derive: 



Gf e f M = (W a ,b + Gft B ) « i?T(lnl0 10 - ln^ B ) (14) 

At physiological temperatures the following region of energy corresponds to the 
frequency range of pulsations (13) 



Gref w (4 - 8) kcal/mole (15) 

Such a region of energies is pertinent to a wide range of biochemical processes. 

The quick jump- way pulsations of macromolecules can cause acoustic shock- 
waves in the solvent and its structure destabilization. Concomitant increase 
in water activity leads to distant interaction between different proteins as well 
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as proteins and cells. Such solvent-mediated phenomena were discovered and 
studied in our laboratory by set of specially elaborated methods (Kaivarainen, 
1985, 1986, 1987; Kaivarainen et al., 1990, Kaivarainen et al., 1993). 

When the [A B] transitions in proteins are related to continuous diffusion 
only, then the G s t values calculated using eq.(3) for the same frequency interval 
(10 4 - 10 7 ) s, is about (3 - 7) kcal/mole. 

The Kramers equation (1940), which has earlier been widely used for de- 
scribing diffusion processes, has the form: 



where A is a constant and 77 - solvent viscosity. 

The pre-exponential factor in our generalized equation (3) con- 
tains not only the viscosity, but also the temperature and the effec- 
tive volume of a particle. It was shown in our experiments that eq. 
(3) describes the dynamic processes, which occur in macromolecules, 
solutions much better than the Kramer's equation (16). 



3. Dynamic model of protein-ligand complexes formation 

According to our model of specific complexes formation the fol- 
lowing order of events is assumed (Fig. 2): 

1. Ligand (L) collides with the active site (AS), formed usually by 
two domains, in its open (b) state: the structure of water cluster in 
AS is being perturbed and water is forced out of AS cavity totally or 
partially; 

2. Transition of AS from the open (b) to the closed (a) state occurs 
due to strong shift of [a b] equilibrium to the left, i.e. to the AS 
domains large scale dynamics; 

3. A process of dynamic adaptation of complex [L+AS] begins, 
accompanied by the directed ligand diffusion in AS cavity due to 
its domains small-scale dynamics and deformation of their tertiary 
structure; 

4. If the protein is oligomeric with few AS, then the above events 
cause changes in the geometry of the central cavity between subunits 
in the open state leading to the destabilization of the large central 
water cluster and the shift of the A ^ B, corresponding to R ^± T equi- 
librium of quaternary structure leftward. Water is partially forced 
out from central cavity. 

Due to the feedback mechanism this shift can influence the [a b] 
equilibrium of the remaining free AS and promotes its reaction with 
the next ligand. Every new ligand stimulates this process, promoting 




(16) 
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the positive cooperativity. The negative cooperativity also could be 
resulted from the interaction between central cavity and active sites; 

5. The terminal [protein — ligand] complex is formed as a conse- 
quence of the relaxation process, representing deformation of domains 
and subunits tertiary structure. This stage could be much slower than 
the initial ones [1-3]. As a result of it, the stability of the complex 
grows up. 

Dissociation of specific complex is a set of reverse processes to that 
described above which starts from the [a* — > b] fluctuation of the AS cavity. 

In multidomain proteins like antibodies, which consist of 12 do- 
mains, and in oligomeric proteins, the cooperative properties of H2O 
clusters in the cavities can determine the mechanism of signal trans- 
mission from AS to the remote effector regions and allosteric protein 
properties. 

The stability of a librational water effecton as coherent cluster strongly de- 
pends on its sizes and geometry. This means that very small deformations 
of protein cavity, which violate the [cavity-cluster] complementary condition, 
induce a cooperative shift of [A B] equilibrium leftward. The clusterphilic 
interaction, introduced by us (see section 13.3 of [1] or paper [4]) turns to hy- 
drophobic one due to [Ib/tr] conversion. 

This process can be developed step by step. For example, the reori- 
entation of variable domains, which form the antibodies active site (AS) after 
reaction with the antigen determinant or hapten deforms the next cavity be- 
tween pairs of variable and constant domains forming F b subunits (Fig. 2). The 
leftward shift of [A <*=>■ B] equilibrium of this cavity, in turn, changes the ge- 
ometry of the big central cavity between F a b and F c subunits, perturbing the 
structure of the latter. Therefore, the signal transmission from the AS 
to the effector sites of F c subunits occurs due to the balance shift 
between clusterphilic and hydrophobic interactions. This signal is re- 
sponsible for complement- binding sites activation and triggering the receptors 
function on the lymphocyte membranes, in accordance to our model. 

The leftward shift of [A -o- B] equilibrium in a number of cavities in the 
elongated multidomain proteins can lead to the significant decrease of their 
linear size and dehydration. The mechanism of muscular contraction is probably 
based on such phenomena and clusterphilic interactions (see next section). 

For such a nonlinear system the energy is necessary for reorientation of the 
first couple of domains only. The process then goes on spontaneously with 
decreasing the averaged protein chemical potential. 

The chemical potential of the A- conformcr is usually lower than that of 
B- conformer (Ga < &b) and the relaxation of protein is accompanied by the 
leftward A <^> B equilibrium shift of cavities. 

It is predictable, that hydration of proteins will decrease, when clusterphilic 
[water-cavity] interaction turns to hydrophobic one. 



16 



I 



IV 




ii in 



Fig. 2. The schematic picture of the protein association (Fab 
subunits of antibody with a ligand), which is accompanied by the 
destabilization of water clusters in cavities, according to the dynamic 
model (Kaivarainen, 1985). The dotted line denotes the perturba- 
tion of the tertiary structure of the domains forming the active site. 
Antibodies of IgG type contain usually two such Fab subunit and 
one Fc subunit, conjugated with 2Fab by flexible hinge, forming the 
general Y-like structure. 

Our dynamic model of protein behavior and signal transmission, described 
above, is an alternative to solitonic mechanism of non dissipative signal trans- 
mission in proteins and in other biosystems proposed by Davidov (1973). Prop- 
agation of solitonic wave is a well known nonlinear process in the ordered ho- 
mogeneous mediums. The solitons can originate, when the nonlinear effects 
are compensated by the wave dispersion effects. Dispersion is reflected in fact 
that the longer waves spreads in medium with higher velocity than the shorter 
ones. 

However, biosystems of nonregular, fluctuating structure are not 
the mediums, good for solitons emergency and propagation. 

Our dynamic model takes into account the real multidomain and 
multiglobular structure of a proteins and properties of their hydration 
shell fractions. In contrast to Davidov's solitonic model, the dissipa- 
tion processes like reversible "melting" of water clusters, accompa- 
nied by large-scale dynamics of proteins, are the necessary stages of 
our [hydrophobic ^ clusterphilic] mechanism of signal transmission 
in biosystems. 
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The evolution of the ideas of the protein-ligand complex formation 
proceeded in the following sequence: 

1. "Key-lock" or the rigid conformity between the geometry of an 
active site and that of a ligand (Fisher, 1894); 

2. "Hand-glove" or the so-called principle of induced conformity 
(Koshland, 1962); 

3. At the current stage of complex-formation process understand- 
ing, the crucial role of protein dynamics gets clearer. Our model 
allows us to put forward the "Principle of Stabilized Conformity 
(PSC)" instead that of "induced conformity" in protein-ligand spe- 
cific reaction. 

Principle of Stabilized Conformity (PSC) means that the geome- 
try of the active site (AS), optimal from energetic and stereochem- 
ical conditions, is already existing BEFORE reaction with ligand. 
The optimal geometry of AS is to be the only one selected among 
the number of others and stabilized by ligand, but not induced "de 
nova" . 

For example, the [a b] large-scale pulsations of the active sites 
due to domain fluctuations and stabilization of the closed (a) state 
by ligand are necessary for the initial stages of reaction. Such active 
site pulsations decreases the total activation energy necessary for the 
terminal complex formation as multistage process. 

4. The life-time of quasiparticles and frequencies of their excitation 

The set of formula, describing the dynamic properties of quasipar- 
ticles, introduced in mesoscopic theory was presented at Chapter 4 
of book [1] and paper [2]: 

The frequency of c- Macrotransitons or Macroconvertons excitation, repre- 
senting [dissociation/association] of primary librational effectons - "flickering 
clusters" as a result of interconversions between primary [lb] and [tr] effectons 
is: 



where: Pmc = Pac ■ Pbc is a probability of macroconverton excitation; 
Z is a total partition function (see eq.4.2 of [1, 2]); 
the life-time of macroconverton is: 



F C M = • Pmc/Z 

TMc 



(17) 



TMc = (t oc • Uc) 



1/2 



(18) 
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The cycle-period of (ac) and (be) convertons are determined by the sum of 
life-times of intermediate states of primary translational and librational effec- 
tons: 



T~ac {l~a)tr "i" [Ta)lbi (1Q\ 

n c = (n)tr + (n)ib; 

The life-times of primary and secondary effectons (lb and tr) in a- and b- 
states are the reciprocal values of corresponding state frequencies: 



[T a = l/v a : Ta = l/Ua} tr ,lb\ [n = l/v b ; T b = l/l^]tr,lb (20) 

[(^ a ) and {vb)]tr,ib correspond to eqs. 4.8 and 4.9 of [1, 2]; 

[(va) and (i^)]tr,ib could be calculated using eqs. 4. 16; 4.17 [1, 2]. 

The frequency of (ac) and (be) convertons excitation [lb/tr]: 



— • Pac/Z (21) 



F bc = — - Pbc/Z (22) 

Uc 

where: P ac and Pb c are probabilities of corresponding convertons excitations 
(see eq.4.29a of [1, 2]). 

The frequency of Supereffectons and Superdeformons (biggest fluc- 
tuations) excitation is: 

Fsd = — - , ■ Pf /Z (23) 

It is dependent on cycle-period of Supereffectons: tsd = ta* + tb* + td* 
and probability of Superdeformon activation (Pg ), like the limiting stage 
of this cycle. 

The averaged life-times of Supereffectons in A* and B* state are dependent 
on similar states of translational and librational macroeffectons : 



TA* = [{T A )tr ■ (T A )lb] = [{T a Ta)tr ' (T a Ta)lb} 1/2 (24) 

and that in B state: 

T B * = [(Tfl)tr • [TB)lb] = [{nn)tr ■ (nT-)l b } 1/2 (25) 
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The life-time of Superdeformons excitation is determined by frequency of 
beats between A* and B*states of Supereffectons as: 



T D *=l/\{l/T A *)-{l/T B ,)\ 



(26) 



The frequency of translational and librational macroeffectons A ^ 

B cycle excitations could be defined in a similar way: 



1 



M 



{ta +t b + t d ) 



Pm/Z 



tr,lb 



where: 



and 



{TA)tr,lb = [(t • Ta)tr,lb] ^ 2 



{TB)tr,lb = [{ T b ■ T^trAb} 1 ^ 2 



(27) 

(28) 
(29) 



(T D )tr.ib = l/\(l/r A )-(l/T B )\ trlb (30) 

The frequency of primary translational effectons (a ^b) tr transi- 
tions could be expressed like: 



; ' Z (/',/)„• (31) 



(t q + n + T t 



where: {Pd)tr is a probability of primary translational deformons excitation (eq. 
4.25 of [1, 2]); 

[Ta' : Tb]tr are the life-times of (a) and (b) states of primary translational ef- 
fectons (eq. 20). 

The frequency of primary librational effectons as (a ^ b)ib cycles 
excitations is: 



Fib = 71 , ^ , • (Pd)ib (32) 

(T a +U+ T t )lb 

where: (Pd)z&is a probability of primary librational deformons excitation; r a and 
Tb are the life-times of (a) and (b) states of primary librational effectons defined 
as (20). 

The life-time of primary transitons (tr and lb) as a result of quantum beats 
between (a) and (b) states of primary effectons could be introduced as: 
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[r t = \l/r a -l/T b \ 



tr,lb 



(33) 



For the case of (a <^> &) 1 ' 2 * 3 transitions of primary and secondary effectons 
(tr and lb), their life-times in (a) and (b) states are the reciprocal value of 
corresponding frequencies: [r a = \jv a and = l/vbll^w- These parameters 
and the resulting ones could be calculated from eqs.(2.27; 2.28 of [1]) for primary 
effectons and (2.54; 2.55 of [1]) for secondary ones. 

The results of calculations, using eqs. (31, 32) for frequency of excitations 
of primary tr and lb effectons are plotted on Fig. 3a,b. 

The frequencies of Macroconvertons and Superdeformons were calculated 
using eqs. (17 and 23). 



[v a , b l tr> 10 10 (1/s) 



[v,/bk, 10 s (1/s) 
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Fig. 3. (a) - Frequency of primary [tr] effectons excitations, 
calculated from eq.(31); 

(b) - Frequency of primary [lb] effectons excitations, calculated 
from eq.(32); 

(c) - Frequency of [Ib/tr] Macroconvertons (flickering clusters) 
excitations, calculated from eq.(17); 

(d) - Frequency of Superdeformons excitations, calculated from 
eq.(28). 

At the temperature interval (0-100)°C the frequencies of translational and 
librational macrodeformons (tr and lb) are in the interval of (1.3-2.8)-10 9 s _1 and 
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(0.2-13)T0 6 s 1 correspondingly. The frequencies of (ac) and (be) convertons 
could be denned also using our software and formulae, presented at the end of 
Chapter 4 of [1, 2]. 

The frequency of primary translational effectons [a b]t r excitations at 
20°C, calculated from eq.(31) is v ~ 7- 10 10 (l/s). It corresponds to electromag- 
netic wave length in water with refraction index (n = 1.33) of: 

A = (cn)/i/ ~ 6mm (34) 

For the other hand, there are a lot of evidence, that irradiation of very different 
biological systems with such coherent electromagnetic field exert great influences 
on their properties (Grundler and Keilman, 1983). 

Between the dynamics /function of proteins, membranes, etc. and 
dynamics of their aqueous environment the strong interrelation is 
existing. 

The frequency of macroconvertons, representing big density fluctuation in 
the volume of primary librational effecton at 37C is about 10 7 (l/s) (Fig 3c). 

The frequency of librational macrodeformons at the same tem- 
perature is about 10 6 s _1 ,i.e. coincides with frequency of large-scale 
protein cavities pulsations between open and closed to water states 
(see Fig. 2). This confirm our hypothesis that the clusterphilic in- 
teraction is responsible for stabilization of the proteins cavities open 
state and that transition from the open state to the closed one is 
induced by coherent water cluster dissociation. 

The frequency of Superdeformons excitation (Fig. 3d) is much lower: 

v s ~ (10 4 - 10 5 ) s- 1 (35) 

Superdeformons are responsible for cavitational fluctuations in liq- 
uids and origination of defects in solids. Dissociation of oligomeric 
proteins, like hemoglobin or disassembly (peptization)of actin and 
microtubules could be also related with such big fluctuations. 

5. Mesoscopic mechanism of enzyme catalysis 

The mechanism of enzyme catalysis is one of the most intriguing 
and unresolved yet problems of molecular biology. It becomes clear, 
that it is interrelated not only with a spatial, but as well with hier- 
archical complicated dynamic properties of proteins (see book: "The 
Fluctuating Enzyme" , Ed. by G.R.Welch, 1986). 

The [proteins + solvent] system should be considered as a cooperative one 
with feedback links (Kaivarainen, 1985, 1992). Somogyi and Damjanovich 
(1986) proposed a similar idea that collective excitations of protein structure 
are interrelated with surrounded water molecules oscillations. 



22 



The enzymatic reaction can be represented in accordance with our 
dynamic model as a consequence of the following stages (Kaivarainen, 
1989; [1] ). 

The first stage: 

(J) E b + S ^ E b S (36) 

- the collision of the substrate (S) with the open (b) state of the active site [AS] 
cavity of enzyme (E). 

The frequency of collisions between the enzyme and the substrate, whose 
concentrations are [Ce] and [Cs], respectively, is expressed with the known 
formula (Cantor and Schimmel, 1980): 



^ col - Anr (D E + D S ) ■ N [C E ] ■ [C s ] (37) 

where: r = a E + as is the sum of the enzyme's and substrate's molecular radii; 
N is the Avogadro number; 



kT kT 
D E = - R and D S = - (38) 

- are the diffusion coefficients of the enzyme and substrate; k is the Boltzmann 
constant; T is absolute temperature; rj is a solvent viscosity. 

The probability of collision of (b) state of the active site with substrate is 
proportional to the ratio of the b-state outer cross section area to the whole 
enzyme surface area: 



fc 

P b = J- E -F b (39) 

where: i 7 }, = j^rf is a fraction of time, the active site [AS] spend in the 
open (b) - state. 

So, the frequency of collision between the substrate and (b) state of the 
active site (AS) with account for (40), meaning the first stage of reaction is: 



v b col = v col -Pb^kj (40) 

The second stage of enzymatic reaction is a formation of the primary enzyme- 
substrate complex: 

(//) E b S^ [E a *S} {1) (41) 
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It corresponds to transition of the active site cavity from the open (6) state to 
the closed (a) one and stabilization the latter state by a ligand. 

The rate constant of the [b — ► a] transitions is derived with the Stokes- 
Einstein and Eyring -Polany generalized equation (3): 

where: n is the solvent viscosity; V is the effective volume of the enzyme 
domain, whose diffusional reorientation accompanies the (b — » a) transition of 
the active site [AS]. 

The leftward shift of the [a 44> b] equilibrium between two states of the active 
site is concomitant with this stage of the reaction. It reflects the principle 
of stabilized conformity, related to AS domains movements that we have put 
forward in the previous section. 

The third stage: 



(7/7) [E a S] <s> [E a 'S*] (43) 

represents the formation of the secondary specific complex. This process is 
related to directed ligand diffusion into the active site cavity and the dynamic 
adaptation of its geometry to the geometry of the active site. 

Here the Principle of Stabilized microscopic Conformity is real- 
ized, when the AS change its geometry from (a) to (a*) without do- 
mains reorientation. The rate constant of this stage is determined by the rate 
constant of substrate diffusion in the closed (a) state of the active site cavity. It 
is also expressed by generalized kinetic equation (42), but with other values 
of variables: 



*J- = iexp (-§)=*.„ (44) 

where: 

r s = (v s /k)rf»/T (45) 

is the correlation time of substrate of volume (v s ) in the (a) state of the active 
site; rf n is the internal effective viscosity; G a s is the activation energy of thermal 
fluctuations of groups, representing small-scale dynamics (SS), which determine 
the directed diffusion of a substrate in the active site [AS] closed cavity. 

The directed character of ligand diffusion in AS can be determined by the 
relaxation of a protein structure, due to perturbation of AS domains by ligand. 
The relaxation changes were observed in many reactions of specific protein- 
ligand complexes formation (Kaivarainen, 1985, 1989). 
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The complex formation [pair of domains forming the AS + substrate], fol- 
lowed by these domain immobilization can be considered as an emergency of a 
new enlarged protein primary effecton from two smaller ones, corresponding to 
less independent AS domains or their compact " nodes" . 

We assume that at this important stage, the waves B of the attack- 
ing catalytic atoms (A B ) and the attacked substrate atoms (A B ) start 
to overlap and interfere in such a way that conditions for quantum- 
mechanical beats between them become possible. 

Let us consider these conditions in more detail. 

According to classical statistics, every degree of freedom gets the energy, 
which is equal to kT/2. This condition corresponds to harmonic approximation 
when the mean potential and kinetic energies of particles are the same: 

V = T k = mv 2 /2 w kT/2 (46) 

The corresponding to such ideal case the de Broglie wave (wave B) length is 
equal to: 



3 mv (mkT)V2 1 ( > 

For such condition the wave B length of proton at room temperature is nearly 
2.5A, for a carbon atom it is about three times smaller and for oxygen - four 
times as small. 

In the latter two cases, the wave B lengths are comparable and even less 
than the sizes of the atoms itself. So, their waves can not overlap and, beats 
between them are not possible. 

However, in real condensed systems with quantum properties, in- 
cluding the active sites of enzymes, the harmonic approximation is 
not valid because (Tk/V) <C 1 (see Fig. 5 of [1]). Consequently, the 
kinetic energy of atoms of AS: 



T k < (l/2)fcT and A B > h/(mkT) 1/2 (47a) 

It must be taken into account that librations, in a general case 
are presented by rotational-translational unharmonic oscillations of 
atoms and molecules, but not by their rotational motions only (Coffey 
et al., 1984). 

The length of waves B of atoms caused by a small translational 
component of most probable impulse, related to librations is bigger 
than that related to pure translations: 

[kb = (h/Pl b )} > [Xt r = {h/Ptr)] (47b) 
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Even in pure water the linear sizes of primary librational effectons 
are several times bigger than that of translational effectons and the 
size of one H 2 molecule (see Fig. 7 of [1] or Fig. 4 of [2]). 

In composition of the active site rigid core the librational waves B 
of atoms can significantly exceed the sizes of the atoms themselves. 
In this case their superposition, leading to quantum beats of waves 
B in the [active site - substrate] complex, accelerating the enzyme 
reaction is quite possible. 

In accordance to our model, periodic energy exchange resulting 
from such beats occurs between waves B of the substrate and active 
site atoms. 

The reaction [S* — > P], accelerating by these quantum beats, is the 
next 4th stage of enzymatic process - the chemical transformation of 
a substrate into product: 



(IV) [E a S*} -> [E a P] (48) 

The angular wave B frequency of the attacked atoms of substrate with mass ms 
and the amplitude As can be expressed by eq.(2.20): 



u s = h/2m s A 2 s (49) 

The wave B frequency of the attacking catalytic atom (or a group of atoms) 
in the active site is equal to: 



U3 



cat 



2m c A c 



(50) 



The frequency of quantum beats which appear between waves B of catalytic and 
substrate atoms is: 



uj* = u cat -lu s =h[ — — - - — A ) (51) 
\2m c A c 2m s AgJ 



The corresponding energy of beats: 



E* = E cat - E s = hw* (52) 

It is seen from these formulae that the smaller the wave B mass of the cat- 
alytic atom (m c ) and its amplitude (A c ), the more frequently these beats occur 
at constant parameters of substrate (ms and As). The energy of beats is 
transmitted to the wave B of the attacked substrate atom from the 
catalytic atom, accelerating the reaction. 
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According to our model, the perturbations in the region of the active site are 
accompanied by the appearance of phonons and acoustic deformons in a form of 
small-scale dynamics of protein structure. They provide the signal transmission 
in oligomcric proteins to the central cavity and other active sites leading to 
allosteric effects. 

It is known from the theory of oscillations (Crawford, 1973) that the effect 
of beats is maximal, if the amplitudes of the interacting oscillators are equal: 



A\ « A% (53) 

The [substrate — ► product] transformation can be considered as a result of the 
substrate wave B transition from the main [S] state to excited [P] state. The 
rate constant of such a reaction in the absence of the catalyst (k s ~* p ) can be 
presented by the modified Eyring-Polany formulae, leading from eq.(2.27 of [1]) 
at condition: ex.p(hv p /kT) >> 1 



E p ( E p -E s 
— 6XP ST - 



= exp 



h{v p 



kT 



(54) 



where: E s = hv s and E p = hv p are the main and excited - transitional 
to product energetic states of substrate, correspondingly; v s and v p are the 
substrate wave B frequencies in the main and excited states, respectively. 

If catalyst is present, which acts by the above described mechanism, then 
the energy of the substrate E 5 is increased by the magnitude E cat with the 
quantum beats frequency oj* (51) and gets equal to: 



E Sc = E S + ^cat ( 55 ) 

Substituting E p = hv^ and E Sc = h{v% + v^) in (54), we derive the 
rate constants for the catalytic reaction in the moment of beats (k Sc ^ p ). This 
corresponds to the 4th stage of enzymatic reaction: 



(IV): 



k Sc^P _ ,.P 



= v exp 



h(v p - v s - v c&t ) 



kT 



= fciv 



(56) 



where: v p , v s and v c are the most probable B wave frequencies of the transition 
[S — ► P] state, of the substrate and of the catalyst atoms, correspondingly. 

Hence, in the presence of the catalyst the coefficient of acceleration 
(q) is equal to: 



9cat 



exp 



kT 



(57) 
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For example, at 



hu°/kT w 10; g cat =2.2 • 10 4 (58) 
At room temperatures this condition corresponds to 

E cat = hv cat ~ 6 kcal/molc. 

The beating acts, followed by transitions of a substrate molecule to activated 
by catalyst excited state (S — > 5c), can be accompanied by the absorption of 
phonons or photons with the frequency oj* (51). In this case, the insolation 
of a [substrate - catalyst] system with ultrasound or electromagnetic field of 
the frequency w* should strongly accelerate the reaction when the resonance 
conditions are satisfied. 

It is possible that the resonance effects of this type can account for the ex- 
perimentally revealed response of various biological systems to electromagnetic 
field radiation with the frequency about 6- 10 10 Hz (Deviatkov et al., 1973). We 
have to point out that just such frequency is close to the frequency of (a=^ b) tr 
transitions excitation of primary translational water effectons. The strict cor- 
relation between the dynamics of water and that of biosystems should 
exist on each hierarchic level of time and space. 

Changes in the volume, geometry and electronic properties of the substrate 
molecule, resulting from its transition to the product, change its interaction 
energy with the active site by the magnitude: 

AE'fp = (Eg* — Ep ) (59) 

It must destabilize the closed state of the active site and increase the probability 
of its reverse [a* — > b*] transition. 

Such transition promotes the last 5th stage of the catalytic cycle 
- the dissociation of the enzyme-product complex: 

(V) [E a "P] <s> [E P P] <&E b + P (60) 

The resulting rate constant of this stage, like stage (II), is described by the 
generalized Stokes-Einstein and Eyring-Polany equation (42), but with different 
activation energy G" t ~* fc valid for the [a* — > b] transition: 

k ^ b = w exp {- 9 ^) = 1/T ^ b (61) 

If the lifetime of the (a*) state is sufficiently long, then the desorption of the 
product can occur irrespective of [a* — > b] transition, but with a longer char- 
acteristic time as a consequence of its diffusion out of the active site's " closed 
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state". The rate constant of this process (fc*) is determined by small-scale dy- 
namics. It practically does not depend on solvent viscosity, but can grows up 
with rising temperature, like at the 3d stage, described by (44): 



where 77^" is active site interior viscosity in the closed (a*) state; vp is the 
effective volume of product molecules; P » = cxp(— G P / 1 RT) is the probability 
of small-scale, functionally important motions necessary for the desorption of 
the product from the (a*) state of the active site; is the free energy of 
activation of such motions. 

Because the processes, described by the eq.(61) and (62), are in- 
dependent, the resulting product desorption rate constant is equal 
to: 



k v = k P ^ b +k* P = 1/t p ~ b + l/Tp (63) 
The characteristic time of this final stage of enzymatic reaction is: 



*V = Vkv = ^ (64) 

This stage is accompanied by the relaxation of the perturbed AS domain and 
remnant protein structure to the initial state. 

After the whole reaction cycle is completed the enzyme gets ready for the 
next cycle. The number of cycles (catalytic acts) in the majority of enzymes is 
within the limits of (10 2 — 10 4 ) s _1 . It means that the [a ^ b) pulsations of the 
active site cavities must occur with higher frequency as far it is only one of the 
five stages of enzymatic reaction cycle. 

In experiments, where various sucrose concentrations were used at constant 
temperature, the dependence of enzymatic catalysis rate on solvent viscosity 
{T/rj) was demonstrated (Gavish and Weber, 1979). The amendment for chang- 
ing the dielectric penetrability of the solvent by sucrose was taken into account. 
There are reasons to consider stages (II) and/or (V) in the model described 
above as the limiting ones of enzyme catalysis. According to eqs.(43) and (61), 
these stages depend on {T/rj), indeed. 

The resulting rate constant of the enzyme reaction could be ex- 
pressed as the reciprocal sum of life times of all its separate stages 
(I-V): 



/c ros = l/r rcs = 1/(77 + tii + Tin + nv + t v ) (65) 

where 
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The slowest stages of the reaction seem to be stages (II), (V), and stage 
(III) . The latter is dependent only on the small- scale dynamics in the region of 
the active site. 

Sometimes product desorption goes on much more slowly than other stages 
of the enzymatic process, i.e. 

ky < fcni < fcn 

Then the resulting rate of the process (fc r0 s) is represented by its limiting stage 
(eq. 63): 



n^rcs ~ exp i 



Ga - 
st 



RT 



+ k* P = l/r r( 



(66) 



The corresponding period of enzyme turnover: T ros w ry (eg 64). The 
internal medium viscosity (r) m ) in the protein regions, which are far from the 
periphery, is 2-3 orders higher than the viscosity of a water-saline solvent (0.001 
P) under standard conditions: 

(f] in /f]) > 10 3 

Therefore, the changes of sucrose concentration in the limits of 0-40% at con- 
stant temperature can not influence markedly internal small- scale dynamics in 
proteins, its activation energy (G a ) and internal microviscosity (Kaivarainen, 
1989b). This fact was revealed using the spin-label method. It is in accordance 
with viscosity dependencies of tryptophan fluorescence quenching in proteins 
and model systems related to acrylamide diffusion in protein matrix (Eftink 
and Hagaman, 1986). In the examples of parvalbumin and ribonuclease T\ it 
has been shown that the dynamics of internal residues is practically insensitive 
to changing solvent viscosity by glycerol over the range of 0.01 to 1 P. 

It follows from the above data that the moderate changes in solvent viscosity 
(77) at constant temperature do not influence markedly the k* p value in eq.(66). 
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Therefore, the isothermal dependencies of k ros on (T/rf)T with changing 
sucrose or glycerol concentration must represent straight lines with the slope: 



Afc ros k ( G a st b \ , 

t9a -WJvW^v c ^{-n^r) T (67) 

The interception of isotherms at extrapolation to (T/rj — ► 0) yields (62): 



kT ( G a * \ 

(fcrcs) (T/l)) ^ = k *P = eX P ) (68) 

The volume of the Brownian particle (V) in eq.(67) corresponds to the effective 
volume of one of the domains, which reorientation is responsible for (a^ b) 
transitions of the enzyme active site. 

Under conditions when the activation energy G" t ^ b weakly depends on tem- 
perature, it is possible to investigate the temperature dependence of the effective 
volume V, using eq.(67), analyzing a slopes of set of isotherms (67). 

Our model predicts the increasing of V with temperature rising. This reflects 
the dumping of the large-scale dynamics of proteins due to water clusters melting 
and enhancement the Van der Waals interactions between protein domains and 
subunits (Kaivarainen, 1985; 1989b, Kaivarainen et al., 1993). The contribution 
of the small-scale dynamics (k*) to fc ros must grow due to its thermoactivation 
and the decrease in rj ln and G a (eg. 32). 

The diffusion trajectory of ligands, substrates and products of en- 
zyme reactions in "closed" (a) states of active sites is probably deter- 
mined by the spatial gradient of minimum wave B length (maximum 
impulses) values of atoms, forming the active site cavity. 

We suppose, that functionally important motions (FIM), introduced by 
B. Frauenfelder et al., (1985, 1988), are determined by specific geometry of the 
impulse space characterizing the distribution of small-scale dynamics of domains 
in the region of protein's active site. 

The analysis of the impulse distribution in the active site area 
and energy of quantum beats between de Broglie waves of the atoms 
of substrate and active sites, modulated by solvent-dependent large- 
scale dynamics, should lead to complete understanding of the physical 
background of enzyme catalysis. 

6. The mechanism of ATP hydrolysis energy utilization in muscle 

contraction 
and protein polymerization 

A great number of biochemical reactions are endothermic, i.e. they 
need additional thermal energy in contrast to exothermic ones. The 
most universal and common source of this additional energy is a re- 
action of adenosinetriphosphate (ATP) hydrolysis: 
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fcl 

ATP & ADP + P 
fc-i 



(69) 



The reaction products are adenosincdiphosphatc (ADP) and inorganic phos- 
phate (P). 

The equilibrium constant of the reaction depends on the concentration of 
the substrate [ATP] and products [ADP] and [P] like: 



k L= [ADP] ■ [P] 

fc_i [ATP] 1 ' 

The equilibrium constant and temperature determine the reaction free energy 
change: 



AG = -RT \nK = AH- TAS (71) 

where: AH and AS are changes in enthalpy and entropy, respectively. 

Under the real conditions in cell the reaction of ATP hydrolysis is 
highly favorable energetically as is accompanied by strong free energy 
decrease: AG = -(11 13) kcal/M. 

It follows from (71) that AG < 0, when 



TAS* > AH (72) 

and the entropy and enthalpy changes are positive (AS > and AH > 
0). However, the specific molecular mechanism of these changes in 
different biochemical reactions, including muscle contraction, remains 
unclear. 

Acceleration of actin polymerization and tubulin self-assembly to 
the microtubules as a result of the ATP and nucleotide GTP splitting, 
respectively, is still obscure as well. 

Using our model of water-macromolecule interaction [6], we can 
explain these processes by the "melting" of the water clusters - libra- 
tional effectons in cavities between neighboring domains and subunits 
of proteins. This melting is induced by absorption of energy of ATP 
or GTP hydrolysis and represents [Ib/tr] conversion of primary libra- 
tional effectons to translational ones. It leads to the partial dehydra- 
tion and rapprochement of domains and subunits. The concomitant 
transition of interdomain/subunit cavities from the "open" B-state 
to the "closed" A-state should be accompanied by decreasing of lin- 
ear dimensions of a macromolecule. This process is usually reversible 
and responsible for the large-scale dynamics. 
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In the case when disjoining clusterphilic interactions that shift the 
[A B] equilibrium to the right are stronger than Van der Waals 
interactions stabilizing A-state, the expansion of the macromolecule 
can induce a mechanical "pushing" force. 

In accordance to our model, this "swelling driving force" is re- 
sponsible for shifting of myosin "heads" as respect to the actin fila- 
ments and muscle contraction. 

Such FIRST relaxation "swelling working step" is accompanied by 
dissociation of products of ATP hydrolysis from the active sites of myosin heads 
(heavy meromyosin). 

The SECOND stage of reaction, the dissociation of the complex: [myosin 
"head" + actin], is related to the absorption of ATP at the myosin active site. 
At this stage the A B equilibrium between the heavy meromyosin conformers 
is strongly shifted to the right, i.e. to an expanded form of the protein. 

The THIRD stage is represented by the ATP hydrolysis, (ATP -> 
ADP + P) and expelling of (P) from the active site. The concomitant 
local enthalpy and entropy jump leads to the melting of the water 
clusters in the cavities, B — ► A transitions and the contraction of free 
meromyosin heads. 

The energy of the clusterphilic interaction at this stage is accumulated in 
myosin like in a squeezed spring. 

After this 3d stage is over the complex [myosin head + actin] forms 
again. 

We assume here that the interaction between myosin head and 
actin induces the releasing of the product (ADP) from myosin active 
site. It is important to stress that the driving force of "swelling 
working stage" : [A — > B] transition of myosin cavities - is represented 
by our clusterphilic interactions (see Section 13.3 of [1] and paper 
[5])- 

A repetition of such a cycle results in the relative shift of myosin 
filaments with respect to actin ones and finally in muscle contraction. 

The mechanism proposed does not need the hypothesis of Davydov's soliton 
propagation (Davydov, 1984) along a myosin macromolecule. It seems that 
this nondissipative process scarcely takes place in strongly fluctuating biological 
systems. Soliton model does not take into account the real mesoscopic structure 
of macromolecules and their interaction with water as well. 

Polymerization of actin, tubulin and other globular proteins com- 
posing cytoplasmic and extracell filaments due to hydrophobic inter- 
action can be accelerated as a result of their selected dehydration 
due to local temperature jumps in mesoscopic volumes where the 
ATP and GTP hydrolysis takes place. 

The [assembly disassembly] equilibrium is shifted as a result of such 
mesophase transition to the left in the case when [protein — protein] interface 
Van-der- Waals interactions are stronger than a clusterphilic one. The latter is 
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mediated by librational water effecton stabilization in interdomain or intcrsub- 
unit cavities. 

It looks that the clusterphilic interactions play an extremely im- 
portant role on mesoscale in the self-organization and dynamics of 
biological systems. 

7. Water activity as a regulative factor in the intra- and inter-cell 

processes 

Three most important factors can be responsible for the spatial 
processes in living cells: 

1. Self-organization of supramolecular systems in the form of mem- 
branes, oligomeric proteins and filaments. Such processes can be me- 
diated by water-dependent hydrophilic, hydrophobic and introduced 
by us clusterphilic interactions [1, 5]; 

2. Compartmentalization of cell volume by semipermeable lipid- 
bilayer membranes and due to different cell's organelles formation; 

3. Changes in the volumes of different cell compartments by os- 
motic process correlated in space and time. These changes are depen- 
dent on water activity regulated, in turn, mainly by [assembly^ disassembly] 
equilibrium shift of microtubules and actin's filaments. 

Disassembly of filaments leads to water activity decreasing due to increas- 
ing the fraction of vicinal water, representing a developed system of 
enlarged primary librational effectons near the surface of proteins (see 
section 13.5 of [1] and [5]). The thickness of vicinal water layer is about 50 A, 
depending on temperature and mobility of intra-cell components. 

The vicinal water with more ordered and cooperative structure than that of 
bulk water represents the dominant fraction of intra-cell water. A lot of experi- 
mental evidences, pointing to important role of vicinal water in cell physiology 
were presented in reviews of Drost-Hansen and Singleton (1992) and Clegg and 
Drost-Hansen (1991). 

The dynamic equilibrium: { I II Ill } between three stages of macro- 
molecular self-organization, discussed in Section 13.5 of [1] (Table 2) and in 
paper [5], has to play an important role in biosystems: blood, lymph as well as 
inter- and intra-cell media. This equilibrium is dependent on the water activ- 
ity (inorganic ions, pH), temperature, concentration and surface properties of 
macromolecules. 

Large-scale protein dynamics, decreasing the fraction of vicinal water (Kaivarainen, 
1986, Kaivarainen et al., 1990) is dependent on the protein's active site ligand 
state. These factors may play a regulative role in [coagulation peptization] 
and [gel 4=> sol] transitions in the cytoplasm of mobile cells, necessary for their 
migration. 

A lot of spatial cellular processes such as the increase or decrease in the 
length of microtubules or actin filaments are dependent also on their self- 
assembly from corresponding subunits (a, (3 tubulins and actin). 
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The self-assembly of such superpolymers is dependent on the [association 
(A) dissociation (B)] equilibrium constant {Ka^b = Kg^ A ). In turn, this 
constant is dependent on water activity (a_y 2 o), as was shown earlier (eg. 13. 11a 
and 13.12 of [1] and [5]). 

The double helix of actin filaments, responsible for the spatial organization 
and cell's shape dynamics, is composed of the monomers of globular protein - 
actin (MM 42.000). The rate of actin filament polymerization or depolymer- 
ization, responsible for cells shape adaptation to environment, is very high and 
strongly depends on ionic strength (concentration of NaCl, Ca ++ , Mg ++ ). For 
example, the increasing of NaCl concentration and corresponding decreasing of 
dH 2 stimulate the actin polymerization. The same is true of a and (3 tubulin 
polymerization in the form of microtubules. 

The activity of water in cells and cell compartments can be regulated by 
[Na + — K + ] ATP— dependent pumps. Even the equal concentrations of Na + 
and K + decrease water activity an 2 o differently due to their different interaction 
with bulk and, especially with ordered vicinal water (Wiggins 1971, 1973). 

Regulation of pH by proton pumps, incorporated in membranes, also can be 
of great importance for intra-cell an 2 o changing. 

Cell division is strongly correlated with dynamic equilibrium: [assembly -o- 
disassembly] of microtubules of centrioles.. Inhibition of tubulin subunits disso- 
ciation (disassembly) by addition of D 2 0, or stimulation this process by decreas- 
ing temperature or increasing hydrostatic pressure stops cell mitosis - division 
(Alberts et al., 1983). 

The above mentioned factors enable to affect the A B equilibrium of cav- 
ity between a and (3 tubulins, composing microtubules. These factors action 
confirm our hypothesis, that microtubules assembly are mediated by cluster- 
philic interaction (see section 17.5 of [1] and paper [5]). 

The decrease in temperature and increase in intra-microtubules pressure lead 
to the increased dimensions of librational water effectons, clustrons and finally 
this induce disassembly of microtubules. 

Microtubules are responsible for the coordination of intra-cell space organi- 
zation and movements, including chromosome movement at the mitotic cycle, 
coordinated by centrioles. 

The communication between different cells by means of channels can regulate 
the ionic concentration and correspondent an 2 o gradients in the embryo. 

In accordance with our hypothesis, the gradient of water activity, 
regulated by change of vicinal water fraction in different compart- 
ments of cell can play a role of so-called morphogenic factor necessary 
for differentiation of embryo cells. 

8. Water and cancer 

We put forward a hypothesis that unlimited cancer cell division 
is related to partial disassembly of cytoskeleton's actin-like filaments 
due to some genetically controlled mistakes in biosynthesis and in- 
creasing the osmotic diffusion of water into transformed cell. 
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Decreasing of the intra-cell concentration of any types of ions (Na + , K + , H + , 
Mg 2+ etc.), as the result of corresponding ionic pump destruction, incorporated 
in biomembranes, also may lead to disassembly of filaments. 

The shift of equilibrium: [assembly <^ disassembly] of microtubules (MTs) 
and actin filaments to the right increases the amount of intra-cell water, involved 
in hydration shells of protein and decreases water activity. As a consequence of 
concomitant osmotic process, cells tend to swell and acquire a ball-like shape. 
The number of direct contacts between transformed cells decrease and the water 
activity in the intercell space increases also. 

We suppose that certain decline in the external inter-cell water 
activity could be a triggering signal for the inhibition of normal cell 
division. The shape of normal cells under control of cell's filament 
is a specific one, providing good dense intercell contacts with limited 
amount of water, in contrast to transformed cells. 

If this idea is true, the absence of contact inhibition in the case 
of cancer cells is a result of insufficient decreasing of intercell water 
activity due to loose [cell-cell] contacts. 

If our model of cancer emergency is correct, then the problem of 
tumor inhibition is related to the problem of inter - and intra-cell 
water activity regulation by means of chemical and physical factors. 

Another approach for cancer healing we can propose here is the 
IR laser treatment of transformed cells with IR photons frequencies, 
stimulating superdeformons excitation and collective disassembly of 
MTs in composition of centrioles. This will prevent cells division 
and should have a good therapeutic effect. This approach is based 
on assumption that stability of MTs in transformed cells is weaker 
and/or resonant frequency of their superdeformons excitation differs 
from that of normal cells. 
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